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Abstract 

Roberge and Weiss showed that for SU (N) gauge theories, phase transitions occur in the presence 
of quark chemical potential which is imaginary. We show that if such a phase transition is of 
first order, then even with dynamical quarks, 't Hooft loops of arbitrary Z(N) charge are well 
defined at the phase boundary. We outline the computation of the 't Hooft loop at very high 
temperature, where semi-classical methods can be used. To leading order in weak coupling, the 
't Hooft loop satisfies Casimir scaling in the pure glue theory, but not with quarks. Because the 
chemical potential is imaginary, typically the interaction measure is negative on one side of the 
phase transition. Using a matrix model to model the deconfining phase transition, we compute the 
properties of the Roberge- Weiss phase transition for heavy quarks. 
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I. INTRODUCTION 



Understanding the nature of the phase diagram of Quantum Chromodynamics (QCD) is 
one of the outstanding problems in nuclear physics. At zero chemical potential, the theory 
can be studied by numerical simulations on the lattice. Generally this is not possible at 
nonzero chemical potential because of the sign problem [1]. 

One way to trying to understand the theory for nonzero chemical potential is to consider 
a chemical potential which is purely imaginary: then there is no sign problem, and numerical 
simulations are possible. Although an imaginary chemical potential is not directly physi- 
cal, results can be related to those for real chemical potential using Fourier and Laplace 
transformations [2]. 

Roberge and Weiss first considered imaginary chemical potential, 0, and showed that the 
Z(N) dynamics of a SU(N) gauge theory in an essential way [2]. They showed that at high 
temperature, there are first order transitions in the <p plane. As the temperature is lowered, 
these transitions can be of second order or crossover, depending upon the detailed nature 
of deconfinement with dynamical quarks. This has been studied in effective theories [3], in 
holographic models [4, 5], and with numerical simulations on the lattice [6, 7]. 

In this paper we show that Roberge- Weiss phase transitions have an unexpected connec- 
tion with pure glue theories. Consider SU(N) gauge theories without dynamical quarks. 
Running the Wilson loop in the temporal direction gives the Polyakov loop. This is an order 
parameter for Z(N) electric charge, whose expectation value vanishes in the confined phase 
and is nonzero in the deconfined phase. Similarly, the two point function of Polyakov loops 
exhibits an area law in the confined phase. 

The response of the theory to Z(N) magnetic charges is given by the 't Hooft loop [8-12]. 
Its behavior is converse to that of the Wilson (Polyakov) loops, which obeys an area law 
only in the deconfined phase. The string tension for the area law of the 't Hooft loop [9, 10] 
is equal to the order-order interface tension for Z(N) interfaces [13, 14]. (Throughout this 
paper we assume that the 't Hooft loop is purely spatial; with or without quarks, temporal 
't Hooft loops do not develop an area law at any temperature [11].) 

Dynamical quarks carry Z(N) electric charge, and so modify the behavior of both loops. 
The Polyakov loop is no longer a strict order parameter, but is nonzero at any temperature 
T / 0. Similarly, one expects that the 't Hooft loop acts as in the deconfined phase of the 
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pure glue theory, and so exhibits an area law. However, it is not clear how to define the 't 
Hooft loop in a theory with dynamical quarks [9, 10]. 

In this paper we show that at a Roberge-Weiss phase transition that 't Hooft loops can 
be precisely defined for arbitrary Z(N) charge. We stress that 't Hooft loops are only well 
defined at the point where Roberge-Weiss transition(s) occur. Nevertheless, we find it of 
interest that at least at isolated points, that it is possible to do so in a gauge theory with 
dynamical quarks. 

Establishing this result is not difficult. In a gauge theory, a global Z(N) transformation is 
an overall rotation of the quark field by a constant phase. This can be exactly compensated 
by a shift in an imaginary chemical potential for the quarks. What is less obvious is how 
the boundary conditions of the Roberge-Weiss transition are precisely equivalent to those 
for a Z(N) interface, and thus, to the 't Hooft loop [9, 10]. 

We compute the behavior of the 't Hooft loop at asymptotically high temperature, where 
the calculation can easily be done using semi-classical techniques. In the pure glue theory, 
at leading order and up to corrections ~ g 3 , one finds that the interface tension associated 
with 't Hooft loops of different charges satisfies Casimir scaling [15, 16]. We find that even 
to leading order, dynamical quarks do not respect the Casimir scaling found in the pure glue 
theory. 

We also consider the thermodynamics of Roberge-Weiss phase transitions. We find that 
the interaction measure is negative on at least one side of the Roberge-Weiss transition. This 
happens because both the chemical potential and the associated quark number densities are 
imaginary. Thus their contribution to the energy density can be negative, and so unphysical. 

It is also of interest to know how the 't Hooft loop behaves for temperatures which are not 
asymptotic. To study this, we consider heavy quarks, where the deconfining and Roberge- 
Weiss phase transitions remain of first order. We use a matrix model [18? , 19], which 
was previously used to study the transition for heavy quarks [20]. We study the theory in 
the limit where the deconfining transition remains of first order with in the absence of an 
imaginary quark chemical potential. 

In Sec. (II) we discuss the global Z(N) symmetries of the pure gauge theory and how 
such symmetries arise with dynamical quarks. We also discuss the relationship to 't Hooft 
loops. In Sec. (Ill B) consider the thermodynamics of Roberge-Weiss phase transitions. 
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II. Z(N) SYMMETRY AND 'T HOOFT LOOPS 



A. General analysis 

We work in imaginary time r at a temperature T, so r : — > 1/T. Under a gauge 
transformation U(x,t), the gluons and quarks transform as 

A„(x,t) ^-U{x,T)D fl U\x,T) , q(x,r) ->U(x,T)q{x,T) , (1) 

D^ = d^ + igA^. 

A special class of gauge transformations are those which are symmetric up to elements 
of Z(N). Consider the matrix e 2nlk ^ N lN', this has determinant one, and so is an element 
of SU(N). Since it is proportional to the unity matrix, we can consider aperiodic gauge 
transitions 

U(x, 1/T) = e 2nik/N U(x, 0) , (2) 

where k is an integer. Gluons must be periodic in r, but this aperiodicity in the C/'s does 
alter the gluon boundary conditions, cancelling between U and W. 

Quarks are fermions, and so in the absence of a chemical potential, are anti-periodic in r. 
The boundary conditions are altered by the presence of a chemical potential, /i. Consider 
the analytic continuation from a real chemical potential, /i, to one which is imaginary and 
proportional to the temperature, 

fi = 27Ti<f)T. (3) 

The quarks now satisfy 

q(x,l/T) = -e 2 ^g(£,0) . (4) 

Because of the change in the fermion boundary conditions, when 0^0 the theory can 
exhibit unphysical behavior. For example, <fi = \ turns fermions into bosons. We shall see 
an example of this unphysical behavior later, when we find that the interaction measure is 
negative on one side of the Roberge- Weiss transition line. 

Since the Z(N) transformation of Eq. (2), is proportional to the unit matrix, for the 
quarks it just alters them by an overall phase. Consequently, when ^ 0, a Z(N) transfor- 
mation can be absorbed by a shift in <fr, 

0^0+^- (5) 
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The action is invariant under this transformation. By construction, the boundary conditions 
are as well. Thus it is manifestly a symmetry of the theory, and in particular of the partition 
function, 



This is the Roberge- Weiss symmetry [2]. 

This periodicity can also be understood topologically. In the pure glue theory, the gauge 
group is SU(N)/Z(N), and so there is a global symmetry of Z(N). When = 0, the gauge 
group is just SU(N), and this global Z(N) symmetry is lost. However, by introducing an 
imaginary chemical potential, 0, for free quarks we gain an extra global symmetry of U(l). 
The coupling of gauge fields reduces this U(l) symmetry to one of Z(N). 

We thus have global Z(N) symmetries both in the pure glue theory and, when 0^0, 
for dynamical quarks. Thus we might expect that there are non-contractible loops which 
measure windings in these global Z(N) symmetries. In the pure glue theory, such windings 
are measured by the 't Hooft loop. Since the transformation of is so intimately tied to 
the Z(N) symmetry of the pure glue theory, it is perhaps not so surprising that the 't Hooft 
loop continues to measure such windings, even in the presence of dynamical quarks. Where 
in the phase diagram is less obvious, and is addressed next. 

B. Semi-classical analysis 

In this subsection we work in the limit of high temperature, where all computations can 
be done semiclassically. Many of our conclusions are more general, but it is useful having 
explicit expressions to understand the relevant physics. 

We first need to understand how to generate arbitrary Z(N) transformations from the 
SU(N) algebra. To do so, consider the diagonal matrix 



There are N — k elements with entry k/N, and k elements equal to —1 + k/N. Thus has 
zero trace and is an element of the Lie algebra of SU(N). It is analogous to hypercharge 
matrix for three flavors. Since 




(6) 




diag (k...k, -N + k, . . . - N + k) . 



(7) 




(8) 



Yfc generates a Z(N) transformation with strength k. 
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The thermal Wilson line is given by 

( f 1/T \ 
L(a?) = exp ig / A (x,t)cIt J . (9) 

Then the path 

2nT 

A = q Y k , (10) 

9 

takes one from the ordinary perturbative vacuum for q = 0, L = ljy, to the kth Z(N) 
transform thereof, for q — 1, L = e 2mk ^ N lpf. 

At the outset it is also useful to consider the matrix Y^-h, where 

Y N - k =±diag(N-k...N-k,-k...-k) ; (11) 

this has k elements with entry 1 — k/N, and N — k elements with value —k/N. By simply 
reordering the elements, though, it is clear that 

Y N „ k = - Y fc , (12) 

and so 

e 2niY N _ k = e 2ni(-k)/N ^ _ (13) 

That is, Y^r-fc generates a Z{N) transformation with charge N — k, which by the additive 
Z(N) symmetry is equivalent to charge —k. This symmetry is useful, because we expect that 
the interface tension for charge k should be equal to that for charge N — k. The matrix Yjy-fe 
makes this manifest, although our final expressions do not obviously reflect this symmetry. 
Classically there is no potential for the g's, but one is generated at one loop order, equal 

to 

V ij'=l / 

The qi are the elements of gY^. The gluon contribution involves the adjoint covariant 
derivative, and so only the differences of the q^ enter, q^ = qi — qj. 

The qij's arise as the gluon energies divided by 2-kT . Since the potential for the g^'s arises 
for a Matsubara sum, then it is periodic in the g^-, q^ — > q^ + 1. Consequently, the potential 
only involves the absolute value of g^, modulo one. Thus for the gluon term we can always 
require that q : — > 1. 

The term independent of q is the just (minus) the pressure of an ideal gas of gluons. For 
the qij, there are k elements = q(k/N), and N — k elements = g(— 1 + k/N). The potential 
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arises from the (absolute value) of the qi — qj's. For these elements, the |% |'s are either 
zero, or q (assuming q > 0). There are 2k(N — k) such terms which give q, and so the gluon 
potential is 

V p 1(g) - V$(0) = ^L- fc(JV - fc) g 2 (l - q) 2 . (15) 

The potential is identical at q = and 5 = 1, which illustrates the Z(N) symmetry of the 
pure glue theory between the ordinary vacuum, k = 0, and the kth Z(N) transform. 

In all, the potential is proportional to k(N — k). This is known as Casimir scaling. It is 
satisfied in the pure glue theory up to corrections ~ g 3 times the leading order term [15, 16]. 

For a single, massless flavor, quarks contribute to the potential for q as 

N \ 

-^E^-N) 2 • (16) 

8=1 / 

The quark contribution involves the covariant derivative in the fundamental representation, 
so it is the q^ and not q^ — qj, which enter. Like the gluon potential, as the quark potential 
arises from a sum over the Matsubara frequencies, each q$ enters only as the absolute value, 
modulo one. Thus it is necessary to be careful about the range of the Even so, clearly 
Eq. (16) is invariant under two symmetries. As only the absolute value of the qi enter, one is 
qi — > —q,i. The second is apparent from the form of the potential: assuming that 1 > q^ > 0, 
the quark potential is also invariant under qi — > 1 — q^ 
For the q^, there are N — k elements 

and k elements 

The qi are the Euclidean energies divided by 2nT. Thus the first factor of \ is the nT 
which arises because quarks are fermions, with boundary conditions which are anti-periodic 
in imaginary time. The second term, —0, is the contribution of the imaginary chemical 
potential. Lastly, the terms ~ q arise from the background field. 

While we explicitly compute the potential shortly, most aspects of the physics can be 
understood without going into such details. 

First consider vanishing chemical potential, = 0. As discussed above, for the gluons 
the potential is periodic in q, with q = 1 degenerate with q = 0, Eq. (15). This reflects the 
Z(N) symmetry of the pure glue theory. 
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This is no longer true with dynamical quarks. Then g^'s involve factors of kq/N; for 
q — 1, this is k/N, and generates a nontrivial potential. Computation shows the quark 
potential has a higher value when q — 1 than for q = 0. This occurs because quarks in 
the fundamental representation do not respect the Z(N) global symmetry of the pure glue 
theory. 

Looking just at the g^'s shows how the Roberge-Weiss symmetry of Eq. (6) works. Com- 
pare the ordinary perturbative vacuum, where = q = 0, to a state where = k/N, and 
q = 1. When <f> — q = 0, all g^'s equal |. When = fc/iV, and q = 1, all of the g^'s are either 
| or — |. Since only the absolute value enters, both are equivalent to |. 

Thus we see that a state with = k/N is equivalent to the perturbative vacuum, if we 
shift the background field by k/N. As we increase from 0, there is no background field, 
but at some point, it will shift to that with k/N. The natural point for this to happen is 
halfway in between, when = k/(2N). 

Let us consider the g^'s at this midway point. For = k/(2N) and q = 0, all N elements 
are equal to 

• = sfi-^- (18) 



2 V N, 

When = k/(2N) and q — 1, there are N — k elements equal to 



* = 4 (i - £ ) • 



and /c elements equal to 

~2 V~ 

Even without computation one can see that the values in Eqs. (19) and (20) are equal to 
those in in Eq. (18). The q^s in Eq. (19) are related to those in in Eq. (18) by — > 1 — qf, 
the Qi's in Eq. (20) are equivalent to those in Eq. (18) under qi —q^ 
Explicitly, when = k/(2N), the quark potential is given by 




i 

Q 

2 H 



(21) 



The values of the potential are clearly equal when q = and q — 1; note, however, that 
the second term does involve the absolute value of ^ — g. Also, this potential is not simply 
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proportional to k(N — k), and so does not respect Casimir scaling. As argued above, the 
potential will respect a transformation under k — >■ N — k, since then enters, instead of 
Y k . 

We can use these results to discuss 't Hooft loops. In the theory without dynamical 
quarks, Kovner, Korthals-Altes, and Stephanov [9] showed that the 't Hooft loop has a 
simple physical interpretation. Consider a box which is long in one spatial direction, say 
that in the ^-direction. Put a 't Hooft loop of fcth Z(N) charge around the boundary of the 
box at one end of the box, at z = L. Then consider boundary conditions which are q = at 
z = 0, and q — 1 at z — L. These boundary conditions are identical to that for a order-order 
Z(N) interface of charge k [13-16]. The 't Hooft loop at z = L then forces q to jump from 
1 back to 0, so that in all one has periodic boundary conditions. Neglecting this singularity 
shows that the interface tension for the order-order interface is equal to that for the 't Hooft 
loop. 

Our analysis above shows that across the Roberge- Weiss transition point, for (f> = k/(2N), 
that one has identically the same boundary conditions. Across this point, the vacuum jumps 
from q = on the left hand side, to q = 1 on the right hand side. Consequently, at the 
Roberge- Weiss transition point, we can define the 't Hooft loop in precisely the same manner 
as in the pure gauge theory. 

There are some qualifications which apply. In the pure gauge theory the 't Hooft loop 
depends only upon the area of the loop, and not upon its shape. With dynamical quarks, 
the value of the 't Hooft loop depends upon its shape. By construction, we consider a loop 
which is flat. Because quarks carry electric Z(N) charge, 't Hooft loops with a different 
shape have different area terms. 

We stress that the 't Hooft loop can only be defined at the Roberge- Weiss transition 
point. For example, in the ordinary vacuum, = 0, the states with q = and q = 1 are 
not degenerate. Typically, the state with q = 1 is not even extremal, and there is no way in 
which to define the 't Hooft loop [9, 10]. Depending upon the matter content, it is possible 
that the state with q — 1 is metastable [21]. In this case, one can define a 't Hooft loop, but 
the interface tension will have an imaginary part, reflecting this metastability. A necessary 
condition for the Z(N) loop to be non-contractible is if the Z(N) transformed states are 
absolutely degenerate, as across a Roberge- Weiss transition. 

In the pure gauge theory, the 't Hooft loop measures a non-contractible loop for the 
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global Z(N) symmetry. As we argued in Sec. (HA), there is also a global Z(N) symmetry 
if the quarks have an imaginary chemical potential. Since this global Z(N) symmetry of 
dynamical quarks is intimately tied to the Z(N) transformations of the pure glue theory, it 
is natural that the 't Hooft loop continues to measure the winding in Z(N). 

Our results are valid at asymptotically high temperature, where the Roberge-Weiss phase 
transition is manifestly of first order. As one lowers the temperature, the Roberge-Weiss 
transition can become of second order, or even crossover. A single transition line can also 
split into two, as appears to be indicated by numerical simulations on the lattice [6, 7]. 

The 't Hooft loop, and the associated winding in Z(N), is well defined at a Roberge-Weiss 
transition point if it is of first order. It is necessary to be at the Roberge-Weiss transition 
point so that Z(N) transformed vacua are degenerate. The transition needs to be of first 
order for the loop to exhibit an area behavior. For a transition of second order, or crossover, 
the 't Hooft loop does not behave like an area law, but only like the length of the loop. 

The explicit computation of the 't Hooft loop follows standard methods. A Z(N) interface 
has an electric field, which contributes to the action as 



The interface tension is determined as a semiclassical tunneling between q = and q = 1 . In 
the pure glue theory, this involves the sum of the kinetic term in Eq. (22) and the potential 
of Eq. (15). Since each term is proportional to k(N — k), the interface tension is as well, 
and so respect Casimir scaling, at least to leading order in the coupling constant. 

With dynamical quarks, the kinetic term remains as in Eq. (22), but now the potential 
is a sum of the gluon term in Eq. (15) and the quark term in Eq. (21). Since the quark 
potential is not proportional to k(N — k), the interface tension will be a more complicated 
function of k. 

The appearance of Casimir scaling is to some extent an observation about the structure 
of the theory in weak coupling. Computation to corrections ~ g 3 beyond that of leading 
order shows a small violation of Casimir scaling even in the pure glue theory [15, 16]. So 
there is nothing fundamental in Casimir scaling, nor in that it is violated by dynamical 
quarks. The only symmetry principle which must be respected is that for Z(N) periodicity. 
This requires that the interface tension for k is the same for N — k; this is equivalent to — k, 
which as we have argued, is valid. 




(22) 
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C. Minimal path for the 't Hooft loop 

Before discussing how to compute the 't Hooft loop with dynamical quarks, we have to 
be careful in considering the paths which are possible. We have assumed that the path 
for the interface tension is along the direction Y&, but in principle the path could move in 
other directions. For an SU(N) gauge theory, the general path lies in the subspace of all 
commuting generators; this is the Cartan subalgebra, which has N — 1 dimensions. 

For the pure gauge theory one can show that the straight line path is minimal [15, 16]. 
This has deep geometric reasons: the Yj, form the boundary of the Weyl chamber, which 
is the smallest possible region to describe the Cartan subalgebra [18]. The endpoints of Y& 
are the relevant endpoints for the interface tension, so then it is natural that the boundary, 
along Yfc, is the minimal path which connects these two points. 

With dynamical quarks, the structure of the Weyl chamber at the Roberge- Weiss tran- 
sition is more involved. We have not been able to answer this question for an arbitrary 
numbers of colors, and so satisfy ourselves with working out the simplest possible cases, 
working up from N = 2 to N = 4. 

The case of two colors is trivial. There is only one direction, along Yi ~ 03 ~ diag(l, — 1), 
and so the path necessarily lies along Yi. 

The first nontrivial case arises for three colors, where there are two directions in the 
Cartan sub-algebra. For the quarks, the q^s are 



In the standard Gell-Mann notation, the directions are A3 ~ (1,-1,0) and A 8 ~ Yi ~ 
(1, 1, —2), with associated coordinates q 3 and g§. For three colors, there is only one interface 
tension, as k = —2 is equivalent to k — 1. The endpoints of the interface are given by g§ = 
and 1, with q% = 0. The straight line path is along Yi, with one transverse direction, along 



The qi's of Eq. (23) are for quarks, but we can use them for gluons, since only the 
differences of the q^s enter in the gluon potential, through q^ = q% — qj. Numerically we find 
that for the gluon potential of Eq. (14), the straight line path along Yi is minimal. That 
is, for the path where qs 7^ and 03 = 0, for every value of qs the potential is minimal with 
respect to variations in the transverse direction, along q 3 . As noted, this is because Yi is 
the boundary of the Weyl chamber for three colors. 




(23) 



A3. 



11 



Now consider the quark potential at the Roberge- Weiss transition point, = 1/6, using 
the qi's of Eq. (23) in Eq. (16). As in the pure glue theory for three colors, numerically we 
have checked that the minimal path is a straight line along Yi. 

For four colors there are two possible interfaces, k — 1 and k = 2; the associated elements 
of Z(4) are i and —1, respectively. 

For k = 1 we can parametrize the q^s using the usual Cartan generators, 

ft = Q - A U + | diag(l, -1, 0, 0) + | diag(l, 1, -2, 0) + ^ diag(l, 1, 1, -3) . (24) 

The k = 1 interface is from g 15 : to 1, with q 3 = q s = at either end. The straight line 
path is along Yi ~ (1, 1, 1, —3), with only qi§ nonzero. Numerically we checked that the 
straight line path is minimal, both in the pure glue theory and with dynamical quarks at 
the Roberge- Weiss point for k = 1, where = 1/8. 

This exercise also shows that the potential has nontrivial structure. At the Roberge- Weiss 
transition for k = 1, <f> — 1/8, there is a metastable minimum in the potential when g 15 = 0: 
it occurs for q$ = 1, with q^ = 0. It is metastable in all three directions, but tunnels with 
finite lifetime to the usual Z(4) vacua, which are absolutely stable. Such metastable vacua 
are known to arise for these types of potentials [21]. 

Lastly we consider four colors with k = 2. For the q^s we take 

ft = Q - A 14 + I diag(l, -1, 0, 0) + | diag(0, 0, 1, -1) + | diag(l, 1,-1,-1) . (25) 

For the interface with k = 2, q^ = at one end and q^ = 1 at the other, with q2 = q' 2 = at 
both ends. The straight line path is along Y 2 ~ (1,1,— 1,-1). 

To determine stability of a path it is essential to have a parameterization in three indepen- 
dent directions. Our path is along Y2 ~ (1, 1, —1, —1). The diagonal matrix ~ (1,-1,0, 0) 
is a generator for (the diagonal part of) SU(2) in the first two colors, and is obviously trans- 
verse to Y 2 . To determine the remaining direction, one can use brute force: one computes 
the linear combination of (1, 1, —2, 0) and (1, 1, 1, —3) which is transverse to (1, —1, 0, 0) and 
(1,1,-1,-1). The answer, as in Eq. (25), is ~ (0,0,1,-1). This is just the SU(2) type 
generator for the third and fourth and fourth colors of 5*^7(4) . At least after the fact, this 
is obvious. 

By explicit computation, again one finds that the minimal path is a straight line, with 
g 4 7^ and q 2 = q' 2 = 0. This is true both for the pure glue theory, and for the theory with 
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(massless) dynmical quarks at the Roberge- Weiss transition point, <fi = 1/4. We did not find 
metastable minima when k = 2. 

The examples of k = 1 for three colors, and k = 1 and k = 2 for four colors, suggests 
that for SU(N) at the Roberge- Weiss transition point(s), <fi = k/(2N), that the path for 
the associated interface tension is always a straight line along Y^. At present, we can only 
suggest this as a conjecture and have no general proof. If true, surely it is due to the nature 
of the Weyl chamber with dynamical quarks at the Roberge- Weiss transition point (s). 

If the conjecture of a straight line path is true, it is then immediate to compute the 
't Hooft loop. As for the pure glue theory, one has a path in one dimension, along Y&. 
One is left with a tunneling problem in this one dimension, and it is easy to solve this 
tunneling problem by using "energy" conservation for the associated problem in quantum 
mechanics [14]. This involves both the potential and the kinetic term for the gluons, Eq. 
(22). For the pure glue theory, as both terms are ~ k(N — k), and the interface tension 
follows immediately. With dynamical quarks, since the quark contribution to the potential 
is more involved, even with the conservation of energy one is left with a single integral over 
q which needs to be computed numerically. Certainly the interface tension with dynamical 
quarks does not satisfy Casimir scaling. 

III. THERMODYNAMICS OF ROBERGE- WEISS TRANSITIONS 
A. High temperature 

The total pressure p and the entropy density s are a sum, 



where p g and s g are the gluon contributions, and p q and s q the quark contributions. 

On the left side of the Roberge- Weiss transition, (ft = k/(2N), the gluon contribution is 



P = p g +p q , 



s = 



Sg + S, 



(26) 




(27) 
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while the quarks contribute 



n 2 NN { T A 



Pq 



7 n fh 



H rp 



\m.(n q 



2-irNNfT 3 



2 J k 



(28) 



Here kq/N — is defined between —1/2 and 1/2. The quark number density is imaginary 
because the chemical potential is. 

To define the energy density, we take the standard expression for a real chemical potential, 
and assume it remains valid for an imaginary chemical potential: 



e = — p + sT + fin q = —p + sT — 2nT Im(n 9 ). 



(29) 



Now consider how thermodynamic functions change on either side of the Roberge- Weiss 
transition, <fi = k/(2N). On the left side of the transition, q = 0, while on the right, 
q = 1, so kq/N — = ^fk/(2N) changes sign across the transition. The pressure and the 
entropy density are even in kq/N — (f) and do not change. The quark number density is odd 
in kq/N — (j) and so changes sign. This change in sign for the (imaginary) quark number 
density is the only reason why the energy density changes at the Roberge- Weiss point and 
makes the transition of first order. 

If we assume that <fi and T are both fixed, at asymptotically high temperature the con- 
tribute of the quark number density to the internal energy density is: 



lm(n q 



e p 
— r 3 — r — 2tv(l. 

rp4 rp^ r rpz 



(30) 



Thus at high T the interaction measure is due entirely to the contribution from the quark 
number density, 



A 



e — 3p 

y4 



-2tt< 



lm(n q 



T 3 ' 



(31) 



The interaction measure is nonzero even at high T because we assume that the (imaginary) 
chemical potential is proportional to temperature, [i = 2niT(p. As discussed, it also flips 
sign across the transition. This holds not only at high T, but persists down to temperatures 
close to the transition temperature, as we see in the model calculations which follow. 
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B. Non-perturbative models 



The perturbative potential for the g's is given by Eq. (14). It involves the function 
V 2 (x) = x 2 (l - |x|) 2 , and g y = (q { - qj) mod i- 

The minimum of this potential is always the usual perturbative vacuum, or a Z(N) 
transform thereof. To model the transition to deconfinement, we can add, by hand, non- 
perturbative terms [18? -20]. The involves one new function, V\{x) = \x\(l — \x\), where 
like V^x), this function is defined to be periodic in x, modulo one. 

For all values of N, a successful fit was obtained with the form 

V7 = ~ T 2 Tl £ (-| Kfa) - c 2 V 2 ( % ) + c 3 ) (32) 

The parameters c\ and c 2 are assumed to be independent of temperature, while the temper- 
ature dependence of c 3 is simply 

C 3(T) = C3 (oo) + ( C 3(r d )-C3(oo))^ . (33) 

Driven by the lattice data, most terms in the non-perturbative potential are ~ T 2 . The 
temperature dependence of c 3 also incorporates a MIT bag constant ~ c 3 (T d ) — c 3 (oo). 

The parameters of the model are fit by the transition in the pure glue theory, where we 
assume = 270 MeV. At the outset, the model involves four parameters, ci, c 2 , c 3 (Trf) and 
03(00). One then requires two conditions. First, that the transition occurs at T c . Second, 
that the pressure vanishes at Tj. The first is a reasonable condition at any N. The second is 
an approximation, modeling that the pressure is suppressed by 1/iV 2 in the confined phase, 
relative to that in the deconfined phase. This leaves two parameters, which are fit by the 
value of the latent heat, and the fall off of the pressure with temperature at asymptotically 
high T. 

To include quarks, we follow Ref. [20] and simply add the one loop term for quarks in 
the background field. This is then the only way in which the (imaginary) chemical potential 
enters. For constant quark mass, this contribution is 

V q = 2Ttr Lc J ^[ln(l + e-^- l27rT ^») + l n (l + e -^ +i27rT ^»)] , (34) 



where trf c represents the trace over flavor and color and Ef(p) = yp 2 + m 2 . We assume 
that the up and down quarks are isospin symmetric. We thus adopt the notation that the 
light (up and the down) quark masses are mi, and the strange quark mass is m s . 
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For light quark masses there is a back from to the gluon potential for chiral symmetry 
breaking or restoration, but we can neglect this here. Such back reaction may be discussed 
by using the functional renormalization group or the gluon and ghost potentials from the 
Landau gauge gluon and ghost propagators [22]. When we consider the small quark mass 
region, we should also consider the meson and also the baryon contributions, but inclusion 
of these effects is much more involved, and will be treated later. 

We shall work in the limit of heavy quark masses, where the light and strange quark 
masses are constant, and such effects can be neglected. Therefore, as a first step to construct 
the reliable model of QCD, we investigate the upper part of the Columbia plot at the 
Roberge- Weiss endpoint. 

For later convenience, we introduce the modified Polyakov-loop as 



These quantities are useful because they are periodic under the Roberge- Weiss symmetry. 

Previously it was found the the pressure and other thermodynamical quantities exhibit 
unphysical behavior below Td in the matrix model [20]. This is because the behavior of color 
singlet quantities, such as glueballs, are not included self-consistently. To handle this, we 
modify the potential as 



where q c = (1/3,-1/3,0). This follows a similar modification in Ref. [20] where the in- 
dependent part of the potential was modified in one particular direction in color space. The 
above form is more natural, and suppresses it in a color symmetric manner. 

C. Numerical results 

For a example for the T-dependence of the imaginary part of the modified Polyakov-loop 
(Im \]/), we show the mi = m s case. The right (left) panel of Fig. 1 shows the T-dependence of 
Im\l/ for the matrix model for deconfinement (logarithmic Polyakov-loop effective potential). 

We can see that there is the tri-critical point in the matrix model for deconfinement, 
but does not in the logarithmic Polyakov-loop effective potential down to mi = m s = 1 
GeV. There is the possibility that the logarithmic Polyakov-loop effective potential has the 



(35) 



V(q)^V(q)-V(q c ) + V(q c )W , 



(36) 
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8.2 0.3 0.4 8.2 0.3 0.4 

T [GeV] T [GeV] 

FIG. 1. The T-dependence of Im \E'. Right (left) panel denotes the result of the matrix model 
for deconfinement (logarithmic Polyakov-loop effective potential). The dotted and dashed lines are 
mi = m s = 1 and 2 GeV, respectively. 

tri-critical point in smaller mi and m s region, but we should consider the spontaneous chiral 
symmetry breaking in this case and it is out of our present treatment. 

The Columbia plot at = 1/6 of the matrix model for deconfinement is shown in Fig. 2. 



> 




1 2 3 

m na [GeVI 



FIG. 2. The Columbia plot at the RW endpoint. 

The tri-critical point is determined that the gap of the imaginary and real part of the 
modified Polyakov-loop vanished at critical T. The first-order region is larger than that 
at zero chemical potential [20]. It is consistent with recent LQCD prediction. At least in 
the recent LQCD data [7, 23] suggests that there is the tri-critical point at certain T with 
two-flavor quarks and thus our matrix model matches this prediction. The critical mass 
divided by critical temperature (m c /T c ) are summarized in Table. I. The difference between 
T c and Td is 1 — 2 % in the case of mi = m s . These values are well reproducesd the recent 
LQCD with the strong coupling and hopping parameter expansions [24]. On the other hand, 
we can not find the critical mass for the logarithmic Polyakov-loop potential at least down 
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mi = oo 


mi = m s 


m s = oo 


Ref. [24] 


5.56(3) 


6.66(3) 


6.25(3) 


This model 


5.0 


6.4 


5.9 



TABLE I. Summary for m c /T c in the recent LQCD simulation with the strong coupling and 
hopping parameter expansions [24] and the matrix model for deconfinement. We show our result 
down to first decimal place. 

to mi = rn s = 1 GeV. It is the reason why the first-order region of the Columbia plot at 
= 1/6 in Ref. [25] is too wide. 

The pressure, internal energy density, entropy density and quark number density with 
m c /T c = 6.4 are shown in Fig. 3. 




5 10 5 10 5 10 5 10 

T'T„ T/T„ T/T„ T/T„ 



FIG. 3. The T/T^-dependence of the pressure, energy, entropy and quark number densities with 
m c /T c = 6.4 where = 0.27 GeV. The solid and dashed lines represents those at <p = 1/6 — e and 
4> = 1/6 + e, respectively. 

To investigate the model dependence, the interaction measure(31) is interesting quantity. 
The actual result is shown in Fig. 4. We can see that A has the two-peak structure at 
= 1/6 + e. Both lines finally reach the same absolute value with opposite sign as explained 
in III A. 
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5 

T'T d 



10 



FIG. 4. The T/T^-dependence of the interaction measure. The definition of the solid and dashed 
lines is same as Fig. 3. 

D. Summary 

In this paper we showed that 't Hooft loops of arbitrary Z(N) charge are well defined 
even with dynamical quarks, at least at the phase boundary of a Roberge- Weiss transition 
which is of first order. To leading order in weak coupling, the 't Hooft loop satisfies Casimir 
scaling in the pure glue theory, but not with quarks. 

For three colors we computed thermodynamic behavior at large quark mass using an 
effective matrix model for deconfmement, and considered the form of the Columbia plot. 
The matrix model well reproduces results obtained by lattice QCD with strong coupling 
and hopping parameter expansions for m c /T c , while models with a logarithmic Polyakov- 
loop potential do not. We computed the interaction measure about the Roberge- Weiss 
transition, <p = 1/6 ± e, and find that it is enhanced on one side versus the other. 

Considering an imaginary chemical potential is clearly useful to discriminate between 
various effective models of deconfmement. The relationship to the 't Hooft loop suggests 
that it probes more fundamental aspects of the dynamics in unforeseen ways. 
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